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^ . Abstract 



We start from an MIT-bag model calculation which provides information 
about the constituent quark distributions in the nucleon. The constituent 
quarks, however, are themselves considered as complex objects whose par- 
tonic substructure is resolved in deep inelastic scattering. This gives rise to 
structure functions of the constituent quarks which, in the unpolarized case, 
are fitted to data at a fixed scale employing three model parameters. Using 
<5^-evolution equations the data are also well described at other scales. For 
the spindependent structure functions g\'^ we additionally have to introduce 
polarization functions for valence and sea quarks which are determined by 
exploiting the x-dependence of the available proton data only. A negatively 
polarized sea in the range x > 0.01 is suggested. We are then capable of 
predicting the shape of the neutron structure function <?" which turns out 
to be in good agreement with experiment. Finally we present an estimate 
for the transversely polarized structure function g2j offering the possibility of 
extracting the twist-3 contribution and rating its importance. 



1 Introduction 

During the last two decades the extensive performance of high energy experiments 
has made an important contribution to our understanding of the nucleon substruc- 
ture. In particular deep inelastic scattering (DIS) of leptons from nucleon targets 
enables the determination of the nucleon structure functions which contain basic 
information about the quark and gluon dynamics in hadronic matter. In order 
to study this dynamics on the theoretical side, one should actually exploit QCD 
as the underlying theory of strong interaction physics which, however, is hardly 
understood in the nonperturbative domain, thus offering no possibility of calculating 
structure functions from first principles. For this reason one is forced to utilize 
phenomenological approaches. But even then we face the problem that all available 
models are formulated in terms of few quark degrees of freedom and therefore tailored 
to reproduce the low energy properties of the nucleon, whereas for the evaluation of 
realistic structure functions the nucleon wave function should reflect the dynamics 
of a large number of pointlike partons resolved in DIS. Accordingly the structure 
functions obtained in such models have to be associated with a low momentum scale 
and one has to think about procedures to bring them in relation to reality, i.e. the 
experimental situation at Qexp ^ ^ GeV^. 

In this work we employ the traditional MIT-bag model [|I| to explicitly calculate 
structure functions in the low-Q^ domain 0, |^, ^. The picture we then have in 
mind is to regard the bag quarks as effective objects, in the subsequent sections 
called constituent quarks, which reveal their complicated substructure in DIS. This 
gives rise to distribution functions characterizing the dynamics of the quark par- 
tons and gluons building up these effective quarks [^ . Later in the discussion these 
distribution functions are denoted by (pvi'Ps and (j)g. While the constituent quark 
distributions obtained in the bag model provide the nonperturbative input of the 
structure functions, the parton distributions inside the constituent quarks are con- 
nected with large Q^ as involved in DIS. In this sense the constituent quarks can 
be viewed as a bridge between the nonperturbative and the perturbative regime, 
leading to the convolution model for the substructure of the nucleon 0. 

The first goal of our investigations is to demonstrate that a satisfactory and 
consistent description of the unpolarized structure functions can be achieved in such 
a convolution model approach. In this framework the Q^-evolution equations can 



be completely expressed in terms of the parton distribution functions </), exposing 
the change of the substructure of the constituent quarks with increasing resolution. 
We are then going to apply the formalism to the longitudinally polarized structure 
functions of the nucleon, i.e. Qi'"', where a careful analysis of the latest data 0, 
1^ is presented. In doing this, we solely exploit the x-dependence of gi and use 
neither the extrapolated integrals over the data nor information from hyperon /5- 
decays. Our special attention in this analysis is given to the polarization of the 
sea quarks which turns out to be essential for the understanding of gi. With the 
sea polarization in hands a prediction for the neutron structure function g^ can be 
made and compared with the available data [^, |I^. Finally we will also present 
an estimate for the structure function (72 which will be experimentally accessible 



via transversely polarized scattering events 111], [T2[ and for the first time offers the 
possibility of extracting higher-twist contributions in DIS [|^ . 

The paper is organized as follows. In Sect. 2 we will start with a brief review of 
the DIS formalism and the definitions of structure functions. Sect. 3 deals with the 
interpretation of structure functions resulting from quark model calculations. This 
leads to the concept of the convolution model, as explained in Sect. 4. Thereupon 
we will successively devote us to the discussion of F2, gi, and g2, corresponding to 
Sects. 5, 6, and 7, respectively. In Sect. 8 we give a conclusion of our results as well 
as an outlook on problems which remain to be solved. 

2 Deep inelastic scattering and structure func- 
tions 

The Feynman diagram which describes the deep inelastic lepton-nucleon scattering 
process i+N —^ i' + X in leading order of the electromagnetic interaction is depicted 
in Fig. 1. For further considerations we go into the target rest frame in which 

p=(M,0), k = {E,k), k' = {E',k') (1) 

holds and write the frequently appearing quantities as 

ly = Pll = E-E', (2) 

Q' = -q\ (3) 



, = ^L = ^ (4) 

Employing the usual Feynman rules to evaluate the inclusive differential cross sec- 

rfV a^ E' 



tion, one finds 114 



While the leptonic tensor L^*^ can be calculated explicitly due to the pointlike nature 
of the lepton, the hadronic tensor W^j^u contains all the relevant information about 
the substructure of the nucleon and therefore represents a much more complicated 
object. For the spin averaged case we get two completely symmetric tensors |l^ 



L's = \ll [u{k\ s')Yuik, s)] [uik', s')Yuik, s)] 



s,s 

1 



^Ti[Yi/i + m)Yi/i' + ni)] , (6) 



'^" AttM 



-Y. J d'^e^'^<{p,a\[UO,M0)f\P,^) , (7) 



Jfj_ being the operator of the hadronic electromagnetic current. The unpolarized 
structure functions Wi 2 are defined by means of the most general symmetric tensor 



respecting current conservation and Lorentz invariance |n| 



ws u/ / 2N / , Qf^QA , W2iu,q^) ( p-q \( p-q \ .^. 

W^u = Wi{u, q ) \^-g^, + -^ j + ^, \^p^ - -^%j [P^ - -^^-j ■ (8) 

At that point one usually switches over to dimensionless structure functions Fi^2 
which are related to Wi^2 in the following way: 

Fi(x,Q2) = Ml^i(z/,g2), (9) 

F2(x,g2) = uW2{iy,q^). (10) 



In the naive quark parton model |T6[ these functions can be expressed in terms of 
quark distribution functions qf of the various flavours i in the nucleon N, 

2nf 

F^{x) = 2xFf (x) = 5^6,2 xgf(x) , (11) 

4 = 1 

where the independence on Q^ is a consequence of the neglect of QCD-corrections 
to the dominant photon-parton scattering process. 

In the case of spindependent lepton-nucleon scattering the information which is 
specific for the polarization effects enters the completely antisymmetric parts of the 



tensors ITl, 113 



L^ 



-Tr 



W. 



flU 



AtcM 



d'^e'^Hp,a\[J,iO,Ju 



le 



ixvap 



q^ 



sPMGi+ [sP 



,p-q 



p^ 



s ■ q 



G, 



(12) 
(13) 
(14) 



M "^ M 

with s^ denoting the spin vector of the nucleon. The last expression again constitutes 
the most general ansatz and thus serves as a definition for the functions 6*1,2 which 
are related to the commonly used structure functions gi^2 by 

gi{x,Q^) = M^uG,{iy,q^), (15) 

^2(x,Q2) = Mu^G2iu,q^). (16) 



In the naive parton model one finds []TB| 



g^ix) ^ g^ix) 



9ii^) 



92ix) 



I 2ny. 



l^el (g/^(x) -grw) = ^Ee? Agf (x), (17) 



i=l 



1 ^"^' /. 



^x) - g,^ {x^ 



-^Ee?Agf(x),(18) 



where ql characterize the probabilities to find in a longitudinally polarized nu- 
cleon a quark with fiavour i and spin alignment parallel (antiparallel) to the nucleon 
spin, while q^ are the correponding probabilities for transversely polarized nu- 
cleons. 

In addition to the fundamental process 7*g — »• q the framework of QCD supplies 
further processes which lead to corrections of the naive parton model. Taking such 
modifications into account, the above introduced quark distributions acquire a Q^- 
dependence. Furthermore the QCD-improved formalism also involves the gluon 
distribution g{x,Q^). Adopting the variable 

'ln(ll) 



2 , 
/5o lln(f^) 






(19) 



with /?o = 11 ~ l^f, and using the non-singlet (NS) and singlet (S) quark distribu- 



tions 



g^'(x,Q^) = 




qHx,Q') - 





i=l 



(20) 
(21) 



the leading order QCD-evolution equations finally can be written as |T7 



1 



^NS/^ ,n2\ _ f "'V TD (^ \ ^NS/„, r^2\ 



12\ 



,(23) 



^g^^(x,g^) = y^n^,J-U^^(l/,g^), (22) 

X ^ ^ ^ ^ 

P being the splitting functions of the various splitting processes. These evolution 
equations enable the determination of F2 at any perturbatively accessible scale Q^, 
provided one knows the function at a reference scale Q^. Whereas on the one hand 
for the longitudinally polarized distribution functions analogous equations hold with 
q and g replaced by Ag and Ag, respectively, a careful analysis of the operator struc- 
ture of g2 fl^ on the other hand reveals that g2 also receives twist-3 contributions 
so that its Q^-evolution remains a yet unsolved problem. 

3 Interpretation of quark model calculations 

So far the employment of the parton model merely resulted in a parametrization 
of the structure functions in terms of quark distributions which contain nonpertur- 
bative information about the hadronic bound state. From the theoretical point of 
view it must be the goal to calculate the structure functions, and with those also the 
quark distributions, using appropriate models of the nucleon. Basically this is pos- 
sible with the help of Eqs. (I^l^dl) and (0),(0) after having performed an operator 
product expansion of the current commutator. Equipped with this formalism the 
final step towards the explicit calculation of structure functions consists in inserting 
a nucleon wave function into the matrix elements of (|^) and (0). In order to get re- 
alistic results, this wave function should be formulated in terms of pointlike partons 
which are resolved in DIS. In practice, however, there is so far no knowledge what- 
soever how to transform the dynamics of the elementary constituents of the nucleon 
into a wave function. This is why one is for the presence forced to utilize well estab- 
lished, phenomenological low energy models of the nucleon which are based on few 
quark degrees of freedom. One of these models is the MIT-bag model |jl[ which is 
formulated relativistically and incorporates confinement in a simple and transparent 
manner, thus being very well suited for the calculation of structure functions. Using 
the well known SU{6) nucleon wave functions, the results for the physical structure 



functions in the bag model read as follows [^ ^ Q 

MR e^ { 



Frix] 



27r s"^ — sin e 



d(3(3 



^MRx-£\ 



T^{e,P)+T^{e,P) 



--{e-MRx)To{e,f3)T,{e,(3) 



+ [x 



—X] 



(24) 



9i{x] 



5 MR 



9 2-K e'^ - sin^ e 



dpp 



MRx-e\ 
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T^ie,P) 



MRx 



P 



l]T^ie,P) 



--{e-MRx)To{e,l3)T,{e,(3) 



+ [x 



—X) 



(25) 



9t{x] 



5 MR 



9 2n e^ 



sin^e 



dpp 



'e-MRx 



MR.x~e 
2 
2 



T^ie,P) 



Ttie,P) 



+ (x 



—X] 



with 



T„(e,/3)= dzz^ jn{ez)jn{(3z) 



(26) 



(27) 





In the case of the neutron we have {F^{x))t^ = | {Fi{x))^ , wheras both (diix))-^ 
and {gT{^))hag vanish identically. The quantities e and MR entering the expressions 
for the structure functions denote the lowest bag frequency and the product mass 
times bag radius, respectively, and are completely fixed by the boundary and stabil- 
ity condition of the bag, giving e = 2.04 and MR = 4£ 0, so that Eqs. (||)-(|2|) 
do not contain any free parameter. As a consequence of broken translational invari- 
ance in the static approximation of the bag model the structure functions receive 
nonvanishing support from the unphysical region x > 1. For this reason modified 



versions of the bag model [pO| , pl| are equipped with a Peierls-Yoccoz projection [Q 
onto momentum eigenstates which weakens but not completely cures the support 
problem. In this work we will disregard such a projection since we consider the latter 
as a minor deficiency. In fact, the numerical contribution from the region a; > 1 is 
very small and the normalizations 



dx{Fl{x] 



bag 



1 

2 



(2^ 



dx(gf{x)) = U-] , (29) 

1 

dx(g'^{x)) = (30) 

\ / bag 

are fulfilled in good approximation. Here {gA/gv)ha.g = 1-09 denotes the ratio of the 
weak coupling constants in the bag model which underestimates the experimental 



one by roughly 12% |19 



The results for the structure functions summarized in the previous paragraph 
were obtained with the help of a wave function which is formulated in terms of three 
bag quarks. The dynamics of these degrees of freedom yields a fair description of the 
low energy properties of the nucleon which are characterized by observables such as 
mass, magnetic moment, weak coupling constants, charge radius, and electromag- 
netic form factors. All these quantities are associated with a very low momentum 
scale Q^ ^ IGeV^. Comparing, however, the structure functions derived from the 
nucleon wave function of the bag model with experiment, no agreement can be seen. 
As already mentioned before, this of course is due to the fact that the structure func- 
tions measured in DIS are associated with very large momentum scales Q^ ^ 1 GeV^ 
implying a high resolution of the nucleon substructure. In order to bring then the 
structure functions of the bag model in relation to reality, we regard the quarks of 
the bag model as effective degrees of freedom. To be precise, we understand these 
effective quarks as clusters which in addition to the valence quarks also consist of sea 
quarks and gluons [Q. The effective quarks defined in this sense are designated as 
constituent quarks in the following. Accordingly the structure functions calculated 
in the framework of the bag model, Eqs. (p^)-(|26|), can be interpreted as linear 
combinations of constituent quark distributions, labeled by the index C. In analogy 
to Eqs. ([TT]), (|T3) and (|T|) we then get 



2X\ ^ ^ Vbag 


\ /bag 


^ c 


(31) 




(gn^)) - 

\ /bag 


- i>;e^Ag^(x), 
^ c 


(32) 




\ /bag 


- lyy.Aq-ix). 

^ c 


(33) 



4 The convolution model 

The concept of the constituent quarks serves as a bridge between the nonperturba- 
tive and perturbative regime. At low Q^-scales the resolution is very poor and the 
hadron appears to be a bound state of three constituent quarks which seem to be 
structureless. In reality, however, they have a complicated substructure which is re- 
solved in DIS. On the formal level this picture can be translated into a corresponding 
diagram p3| depicted in Fig. 2. It describes a nucleon N with momentum p con- 



taining a constituent quark C with momentum k which in turn consists of a quark 
parton of flavour i carrying momentum k. This parton absorbs the virtual photon. 
Fig. 2 is a reflection of an impulse approximation which implies the assumptions 
that final state interactions on the one hand between the fragments of the nucleon 
and the constituent quark and on the other hand between the fragments of the 
constituent quark and the parton are neglected. While the latter are suppresed by 
l/Q"^, as in the parton model, there is, however, no argument for the suppression 
of the former because of the absence of a characteristic mass scale. Nevertheless, in 
the weak coupling limit of the constituent quarks, as realized in the bag model, the 
impulse approximation should be justified very well p3| . 

Taking the momentum dependences of the quark distributions, as given in Fig. 2, 
into consideration, the probability of finding a quark i in the nucleon A^ is determined 
by 

gf (P, 9) = E / d'k q^ip, k) qf{k, q) . (34) 

c '' 

After transition to the usual invariants one finally arrives at 



c i y \y / 



(35) 



This is the fundamental relation of the convolution model which describes the parton 
distributions in the nucleon as a convolution of the constituent quark distributions 
in the nucleon and the parton distributions in the constituent quarks. 

Let us now have a closer look at the parton distributions in the constituent 
quarks. Since the formalism has to include valence as well as sea quarks, it is 
advantageous to decompose the functions gf into two corresponding parts |]2^ , 

gf (z, Q') = M^, Q') S^c + '^T^^. Q') • (36) 



While the sea quark distributions have to be characterized by an isospin label in 
general, the valence quark distribution should be isospin independent and has to 
fulfill the normalization constraint 

1 

j dz(t),{z,Q^) = l. (37) 



This together with the 6ic of Eq. (p6D is a refiection of associating with each con- 
stituent quark just one valence quark of the same fiavour. For the sea quark distri- 
butions it is reasonable to impose isospin symmetry with respect to the constituent 
quarks as well as full symmetry of the strange quarks, yielding finally the most 
general relations 

rJ^'iz.Q') = 0:/^(z,Q^) = 0f^(z,Q^)=0f/^(^,Q^), (38) 

rJ''iz,Q') = 0:/^(^,Q2) = 0f^(z,Q2)=0f/^(^,Q2), (39) 

0:/^(z,Q2) = 0f/^(z,Q2) = 0:/^(z,Q2)=0r(^,Q2). (40) 

5 The unpolarized structure function F2 

Based on the relations of the previous section we are now capable of formulating the 
convolution model expressions for the unpolarized structure function F2. In addition 
to the constituent quark distributions q^ this expression also involves the parton 
distributions 0. The former refiect the dynamics of constituent quarks inside hadrons 
and contain information about the confinement mechanism. In the following this 
nonperturbative part of F2 will be described by the bag model calculation (P^. The 
latter contain information about the parton dynamics in the constituent quarks. 
Since, however, there is no way to calculate these functions at a microscopic 
level, one has to use suitable parametrizations at a reference scale Q^. In order to 
minimize the number of free parameters, we start our considerations by equating 
(HI) and (|3|), i.e. 0^/^(2, Q^) = ^u/d^^^q2^ ^ (j)s{z,Q^) and relating (|g) to (|3|) 
by the 2;-independent suppression factor 1/2 [|25|, i.e. (pl^^ {z.Q"^) = |0^/^(z, Q^). 
Bearing the normalization condition (pTl) in mind, we then choose 



2, _ r(a + 3/2) 1 



{z,Ql) = A-n-z)'-. (42) 



These expressions are motivated by Regge arguments and dimensional counting 
rules. Employing ([TT|), (|3TD, (|35D , (|36D and S'f/(6)-symmetry, we finally get 



Frix,Q') 



dyf.^ 



2x I -[Fr{y) 

y \ /bag 



X 



.y 



Q' 



22 

H X 

3 



2/ V 



X 



,2/ 



,Q^ 



(43) 



The total of three parameters entering 0„ and 0s are fitted to experiment at the 
reference scale Ql = 10 GeV^, giving a = —0.17, 6 = 4, and A = 0.137. Fig. 3 shows 
the result in comparison with the data EH . 

With the help of the QCD-evolution equations ( p2D and (^) we can then deter- 
mine the structure function at any other scale Q^ which lies in the perturbatively 
accessible domain. This evolution also involves the gluon distribution which in the 
convolution model can be written as 



g{x,Q'^ 



dy(^ 



6/^(Ff(x)) J-,Q' 

y \ /bag V^ 



(44) 



Since the total Q^-dependence of F2 is fully contained in the functions 0, it is 
desirable to derive evolution equations for them. Performing the transition to the 
moments 



dz z 



n-l 



A 



± 



( A^^ \ 

n 
n 

\ K' ) 



dz z 



n-l 



lu 


z.Q')\ 


u 


z.Q') 


V 'Pgi. 


z,Q') ) 


( P 


*Qg\^) 


P,- 


*gq[^) 


Pg- 


*qq[^) 


\Pg- 


*gg\z) ) 



(45) 



(46) 



(47) 



and taking advantage of the properties of convolution integrals, we finally get the 
solutions 






C^o^-'M^^iQl), 



(4^ 



5\A--Ai 



/30 " 






+2nfAl^ ( C"^^" - r*^" ) Mt{Ql) - Mt{Ql) . (49) 
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The variable ( has been defined in Eq. (|r^), exphcit expressions for the moments A 



can be found in [|T^]. For the gluon distribution which also cannot be deduced by 
means of a microscopic model yet, we use the parametrization 

1 

z 



,{z,Ql)=B-(l-zr, (50) 



with B = 1.48, c = 2.37. According to (44) this parameter choice corresponds to 



g{x, Ql) = 4.27-(l — x)*^ which respects the momentum sum rule /q dxx{q^{x, Q'^) + 
g{x,Q'^)) = 1 in conjunction with the singlet quark distribution and turns out to 
be a reasonable ansatz. Finally the QCD-parameter A entering C, is settled at 
A = 200 MeV. 

Employing standard techniques p^, we reconstruct the distribution functions 



0t,(z, Q"^) and 0^(2;, Q^) from the evolved moments, Eqs. (^8]) , (^), using twelve and 
three moments, respectively. These distributions tell us in which way the parton 
substructure of the constituent quarks changes with the resolution Q"^. The results 
for typical scales are displayed in Figs. 4 and 5. From these we can easily determine 
the Q^-dependence of the structure function Ff; the results are compared to some 
data in Figs. 6 and 7. 

After having restricted ourselves to the unpolarized structure function of the 
proton so far, some remarks about the neutron structure function are in order now. 
As usual we investigate the difference Ff — -F^ which in the convolution model is 
given by 

^(i^f(z/)) <pI-.q' 

y \ /bag 

1 

y y /bag 



[F^,-F-]{x,Q') = \xrMFl{y)) 4^-,Q' 

3 J y \ /bae \y 

H-,Q']-rM-,Q']hi) 

^ i y ^ ^bag L \y J \y J_ 

Assuming (p^^^ = (j)^^^, as before, the sea quark contribution drops and we are left 
solely with the valence part. Fig. 8 shows F2 — F2 in this case for Q^ = 40 GeV^. 
For large values of x we observe an underestimation of the data which is closely 
connected to the limit F2/F2 — > 2/3 as x ^ 1. This is a well known deficiency 
of the S'f/(6)-wave function we used at the bag level. There are attempts to cure 
this problem [^ , but one must introduce at least two additional phenomenological 
parameters at this point, namely the diquark masses of the intermediate spin singlet 
and spin triplet states. Since we are not interested in the limit a; — > 1 anyhow, we 

11 



do not make modifications in this direction. For small x the data are overestimated 
and we get, due to Eqs. ( ]28|) and (|37|), the value 1/3 for the Gottfried sum rule 
GSR = Jq ^{F^{x, Q"^) —F2{x, Q"^)). An excellent description of the low-x data and 
a reproduction of the experimental GSR p9|, ^ can be achieved according to (|5l| ) 



by setting 0"/^ ^ cj)^/^ . This requires, however, the introduction of two parameters, 
namely h^^^ and b^/^ ^ instead of just one parameter h. Since we have not found 
a microscopic mechanism yet which could give rise to a difference between 0**/^ 
and 0"'^, i.e. an S'f/(2)-breaking at the level of the partons within the constituent 
quarks, we will not study the GSR in this context, but defer it to an extended 
version of the convolution model which in addition to the three constituent quarks 
also takes into account meson cloud effects. 

6 The polarized structure function gi 

The investigations of Sect. 5 revealed that the convolution model approach enables a 
satisfactory and consistent description of the unpolarized structure functions, both 
as a function of x and Q^. In the following we are going to apply the formalism 
to the spindependent structure function gi which is given by the parton model 
representation (|17|). Writing the polarized quark distributions as A^^, = Vv(pv and 
= V.,d,„ with 



^V1 



z,Q')-<PUz,Q' 



v,s{z,Qn= r^ ' — "r'' \ (52) 



employing (|32D and making use of S'f/(6)-symmetry at the constituent quark level, 
we find the convolution model relation 

{ y \ /bag \y J \y J 

This equation also implies the in our opinion reasonable assumption that due to 
the comparatively small mass of the strange-quarks, ml/Q"^ <^ 1, their polarization 
should not be considerably suppressed compared to the polarization of the u- and 
d-quark flavours. For reasons of simplicity we assume the same polarization for all 
three flavours. Nonperturbative effects enter the structure functions again via the 

12 



bag model, Eq. (P^D, the functions (p are the same as before and the spindependent 
information is fully contained in the polarization functions Vv and Vs- 

The particular form of (^) is perfectly appropriate to perform a careful analysis 
of the latest data. First of all we have a look at V^. Imposing the fundamental 



Bjorken sum rule constraint ||31|, modified by the leading order QCD-correction, we 

obtain 

1 

'dx{gf{x,Q')-g^{x,Q' 



1 I 9a 
6 \gv 



- — 1 



asiQ' 



TX 



dy giiy) 



bag 



dzV,{z,Q^)(^,{z,Q'), (54) 



so that 



employing 



ldzV,{z,Q^)Mz,Q^)= ( 



1 



as{Q' 



n 



(55) 
and the condition 



). This relation together with the normalization (| 
\Vv{z,Q'^)\ < 1 essentially determines the shape of the polarization function of the 
valence quarks. Following Carlitz and Kaur [p2| we choose the parametrization 



v.{^,Ql 



l+Pv 



'I 



(56) 
P„ has to be settled 



at a fixed resolution scale Ql = 10 GeV^. In order to fulfill (| 
at 0.035. 

We now proceed to an investigation of the sea polarization. As a starting point 
it is rather instructive to see which shape of gi would be produced under the as- 
sumption of an entirely unpolarized sea, i.e. Vs{z,Ql) = 0. The corresponding 
result is depicted in Fig. 9. It is quite conspicuous that such a scenario leads to 
a sizable overestimation of g^ in the range 0.01 < x < 0.2 where the sea plays an 
important role. This observation suggests a negative sea polarization in this range 
which we want to determine in the following. In order to achieve this goal, we try 
to reproduce the data as well as possible with the help of a parametrization for Vs 
based on an analogy to Eq. (pBf). For this purpose we find the choice 



Vs{z,Ql) = - 



l + Ps 



(l-zf 



-1 



(57) 



to be fairly appropriate. Fig. 10 shows the sea contribution to the polarized structure 
function gi which is, combined with the previously discussed valence contribution. 
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necessary to describe the available data. This analysis can be directly transformed 
into the polarization function of the sea quarks, resulting in Ps ~ 0.25. The reliabil- 
ity of this parametrization is, however, restricted to the range 0.01 < x < 0.3 since 
on the one hand the influence of the sea quarks is absolutely negligible for larger x 
so that in the range x > 0.3 the sea polarization cannot be extracted from the data 
anyway. Moreover, in the region x < 0.01 the SMC-data feature a rapid increase 
and even lie above the valence part of Qi, thus indicating that there may be a sign 
flip in the polarization of the sea quarks which would cause a breakdown of the 
parametrization (57). Due to the huge error bars in the very small-x region such a 
conclusion is, however, premature and one has to admit that there is at present no 
way to exactly determine the behaviour of the polarization function for x < 0.01. 
This is just the crucial range for the evaluation of the integral over gf which sensi- 
tively depends on the extrapolation to x = 0. Therefore we believe that an analysis 
of the spin structure of the nucleon purely based on the first moment of gf is rather 
questionable. 

Taking our parametrization serious even in the limit x — *> 0, we obtain the 
integral /q dx g{{x^ Qq) = 0.140 and for the total quark spin content we find 

S(Q^) ^ Y.JdxAq,,{x,Ql) 

i=l 

= fjdy[g^,{y)^ {jdzV.{z,Ql) 

^^ Vo 

1 

+b j dzVs{z,Ql) 



= 0.60-0.22 = 0.38 (58) 

which is considerably higher than the original results |3^, 0. According to the 




spindependent evolution equations IT^ this decomposition remains unchanged for 
all values of Q^ . This implies the following spin structure of the proton. In the 
low energy regime we have the dynamics of three bag quarks whose relativistic 
motion intrinsically also involves orbital angular momenta so that the quark spin 
contribution reduces from 1 to 0.60 at this level. In a DIS experiment the virtual 
photon resolves the valence and sea quarks inside the constituent quarks. While 
the valence quarks are essentially polarized parallel to the constituent quarks, the 
spin alignment of the sea quarks, at least for x > 0.01, causes a partial shielding 

14 



of the spin carried by the valence quarks. The extrapolation to the very small-x 
regime, i.e. x < 0.01, is, however, absolutely not clear yet so that we do not want 
to rush to conclusions. At any rate we believe that our results for J^ dxgi{x,Ql) 
and S give lower bounds for these quantities. Furthermore we will not specify the 
flavour decomposition of S since this clearly depends on our assumptions for the 
distribution of the strange-quarks as well as their polarization. 

Of course, the formalism of the convolution model can also be applied to the 
polarized structure function of the neutron, g^, as well so that a prediction can 
be made. Proceeding on the assumption of S'[/(6)-symmetry, g^ is identical with 
the sea contribution to gf and thus given by the sea part of ([53|). The data of 
SLAG and SMC |10|] having been taken, however, at an average Q^ of 2 GeV^ 
and 4.6 GeV^, respectively, a Q^-evolution of the sp independent quark distributions 
which we parametrized at Q^ = 10 GeV^ is required. Like in the unpolarized case 
the evolution equations can be expressed in terms of the functions characterizing the 
substructure of the constituent quarks, with the only difference that the moments 
involving the sum of spin up and down polarization, Eqs. (|45|)-(^7D, have to be 
replaced now by the respective ones involving the difference. In order to perform 
the Q^-evolution explicitly, also the gluon polarization function is needed in addition 
to the previously studied polarization functions of valence and sea quarks. For lack 
of any better knowledge we take the polarized gluon distribution function from a 
perturbative QCD calculation which we are going to discuss elsewhere [0]. The 
resulting shape of the structure function g'^{x,Q'^ = 3.3 GeV^) at an average Q^ 
between the SLAG- and SMG-data set is presented in Fig. 11, clearly demonstrating 
that our convolution model prediction for g^ is in accordance with experiment. From 
that one can draw the conclusion that the sea polarization based on the analysis of 
gi is fully consistent with the data of the neutron structure function. 

7 The polarized structure function ^2 

A formal analysis in the framework of the operator product expansion (OPE) reveals 
that in contrast to gi the second spindependent structure function g2 also receives 



twist -3 contributions [|1^, [1^ . Accordingly we write 

g,{x, Q') = g^ix, Q') + g^'{x, Q') . (59) 
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The twist-2 contribution oi g2 can be directly reconstructed from gi, yielding [^ 

g^ix, Q')=l^ gi{y, Q') - giix, Q') , (60) 



whereas the twist-3 contribution is a consequence of the quark-gluon correlations 
in the nucleon and with that also a reflection of confinement, thus representing 
very complicated physics. These twist-3 effects cannot be neglected a priori and 
should occur in a model which incorporates any confinement mechanism. In the bag 
model they are relatively large and stem from the influence of the bag surface which 
simulates confinement in a phenomenological way [|l^. The structure function g2 
involving such complicated physics, the parton model interpretation oi gr = gi + g2, 
see (0), is of course rather questionable so that, being rigorous, not even a basis of 
parametrization is available yet. 

In order to get a first impression of g2, however, we take a pragmatic point of 



view in the following and maintain the parton model relation (18) as an exemplary 
ansatz. According to ( PBD the bag model result ( P^D is then interpreted in terms of 
transversely polarized constituent quark distributions which are not identical to the 
longitudinally polarized quark distributions entering gi, but also receive modifica- 
tions due to quark-gluon correlations, i.e. the simulation of confinement. Starting 
from the constituent quark distributions, one again can make the transition to the 
parton level via the convolution model approach. In order to give an estimate of g2, 
let us now assume that the twist-3 effects are fully parametized by the bag model 
calculation. In this case the substructure of the constituent quarks does not change 
compared to the longitudinal polarization and the parton distribution functions can 
be taken over from the analysis of gi, so that gT is just determined by ( ^3]) after 
replacing {gi{y))hag by {gT{y))hag- From that we finally extract g2. The result is 
presented in Fig. 12 which separately displays the twist-2 and twist-3 contributions. 
The twist-3 part shows a very similar characteristic in x like the absolute value of the 
twist-2 part and both are definitely of the same order of magnitude. For the integral 
over g2 we find the value 0.002 which is very close to the Burkhardt-Cottingham 
sum rule /q dxg2{x) = [38|. There is no evidence so far that this sum rule may be 



violated. 

Concluding the considerations of this section, we would like to emphasize once 
more that our previously given result of g2 is not supposed to be a high precision 
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prediction, but thought as a reasonable estimate based on the convolution model. 
For lack of any feasible scheme for the full Q^-evolution of twist-3 structure functions 
there is also no way to predict g2 at various Q^, provided it is known at a scale Qq- 

8 Conclusions and outlook 

In this work we presented a convolution model approach for the investigation of the 
nucleon structure functions. Such a model is based on the concept of effective degrees 
of freedom describing the low energy properties of the nucleon. In DIS, however, 
their fermionic substructure, characterized by the parton distribution functions 0^, 
and 05, is directly probed. These functions involve three parameters which were 
determined at a reference scale Qo = 10 GeV^ by experimental input. The Q^- 
evolution equations could be formulated at the level of the parton distributions 
0, thus unravelling the substructure of the constituent quarks in dependence on 
the resolution scale. Comparing our results for F2{x, Q"^) at various Q^ with the 
data, we found the picture of the convolution model to be fully consistent with 
experimental findings. In the case of the difference Fj' — -F^ we observed, however, 
obvious discrepancies in the large- as well as in the small-x region which can be 
attributed to the use of S'f/(6)-symmetry at the constituent quark level and isospin 
symmetry of the sea quarks, respectively. 

Subsequent to the investigation of the unpolarized structure functions we applied 
the formalism of the convolution model to the spindependent structure functions 
Qi'"^. In addition to the nonperturbative information of the bag model calculation 
and the distribution functions the expressions for gf'^ also involve polarization 
functions of the partons. An analysis of the x-dependence of the proton data taken 
so far implied that the valence quarks are predominantly polarized parallel to the 
constituent quarks while the sea quarks tend to screen the spin contribution of the 
valence quarks. Having extrapolated our parametrization of the sea polarization 
function, which can be pinned down rather well for 0.01 < x < 0.2, to a; = 0, 
we obtained the quark spin decomposition S = Sj, + S^ = 0.60 — 0.22 = 0.38. We 
would like to emphasize, however, that the small x-behaviour of the sea polarization, 
particularly the one for x < 0.01, and accordingly also the integrated spin content of 
the sea quarks remain essentially undetermined yet. For this reason one definitely 
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cannot exclude the possibility of a totally unpolarized sea, i.e. S^ = 0, at the 
moment. With the polarization functions in hands we could finally make a prediction 
for (?". Comparison with experiment indicated that the sea contribution to the 
proton structure function is absolutely consistent with the full neutron structure 
function. From that one can also conclude that the proton and neutron data together 
are perfectly compatible with the fundamental Bjorken sum rule which is fulfilled 
in our model per construction. 

In the last section of our work we presented an estimate for the transversely 
polarized structure function g2. In doing this we took the naive parton model in- 
terpretation as a reasonable guideline and assumed that the twist-3 effects are fully 
parametrized by the confinement mechanism of the bag model. The resulting es- 
timate for g2 in the convolution model suggests that twist-3 contributions may be 
relatively important even for values of Q^ of the order of 10 GeV^. 

As became clear in the course of our discussions there are some points which 
call for further investigations. The first problem is associated with the small-x re- 
gion of F2 — F2 and consequently connected with the deviation of the GSR from 
1/3. Here we aim at a microscopic mechanism which explains the preference of 
d-quarks over u-quarks in the proton sea. This requires an extension of the convo- 
lution model which in addition to the three constituent quarks also includes a meson 
cloud. Another challenge consists in achieving a microscopic understanding of the 
sea polarization as given by our analysis. The spin problem would then finally loose 
its magic. Especially this issue is intensively worked on [^. 
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Figure captions 

Fig. 1. Deep inelastic lepton-nucleon scattering in the one photon exchange approx- 
imation. 

Fig. 2. DIS in the convolution model for the nucleon substructure. 

Fig. 3. The unpolarized structure function -^^(x, Qq = 10 GeV^) in the convolution 



model. The data are taken from ^B\ 
Fig. 4. Q^-evolution of (pv^z^Q"^). 
Fig. 5. Q^-evolution of 0^(2;, Q^). 



Fig. 6. Ff evolved to Q^ = 5 GeV . The data are taken from ||2^ 



Fig. 7. Ff evolved to Q^ = 50 GeV . The data are taken from ||26 



Fig. 8. [Ff — Ff](x, Q^ = 40 GeV^) under the assumption of an isospin symmetric 
sea. The data are taken from 



Fig. 9. gi{x, Qq = 10 GeV^) in the convolution model under the assumption of an 
identically vanishing sea polarization. The data are taken from [^, ^ Q. 

Fig. 10. Valence and sea contribution to the structure function gi{x, Ql = 10 GeV^). 
The data are taken from f^, |^, |33[ . 



Fig. 11. The neutron structure function g'^{x,Q'^) at an average Q"^ = 3.3 GeV^ in 



comparison with the data M, 10 



Fig. 12. An estimate for g2{x, Qq = 10 GeV^) in the convolution model. 
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